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Abstract 

We establish a nonlinear instability of the Euler-Poisson system for polytropic gases 
whose adiabatic exponents take value in6/5<7<4/3 around the Lane-Emden equilib- 
rium star configurations. 



1 Introduction 

One of the simplest fundamental hydrodynamical models to describe the motion of self- 
gravitating Newtonian inviscid gaseous stars is the compressible Euler-Poisson system: 

cV + V-(pu) = 
p(d t u + u-Vu) + Vp = -pV$ (1.1) 
A<5 = iirp 

where (t, x) £ M + x R 3 and p , u and p denote respectively the density, velocity, and pressure 
of gas. $ is the gravitational potential and it is related to the gas through the Poisson 
equation. We are interested in the polytropic gases (polytropes) for which the pressure is 
assumed to be a power function of the density. More specifically, we consider the following 
equation of state: 

p = Kpi (1.2) 

where K is an entropy constant and 7 > 1 is the adiabatic gas exponent. The values of 7 
distinguish the property of stars. For instance, 7 = 5/3 is used for a monatomic gas, 7 = 7/5 
for a diatomic gas. The smaller 7 represents the heavier molecule of the gas. 

The Euler-Poisson system (jl.ip under the spherically symmetric motion - pit, x) = p(t, r), 
u(t,x) = u(t,r)^ where r = \x\ — reads as follows 

1 / 2 \ 
Pt + — {r pu) r = 0, 

(1.3) 

2. 
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put + puu r + p r -\ 2~ / ps as = 0. 
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Stationary solutions (po(r), uq = 0) of (jl.3p satisfy the following ordinary differential equation 



f + r P s>ds = (1.4) 
dr r z J 

which can be transformed into the famous Lane-Emden equation, and its non-negative so- 
lutions can be characterized according to 7 as follows [2 [12]: Letting M(p) = f A-Ks 2 p{s)ds 
be the total mass of a star, if 7 > 6/5 and any M > 0, then there exists at least one com- 
pactly supported solution p such that M(p) = M. For 7 > 4/3, every solution is compactly 
supported and unique. If 7 = 6/5 and any M > 0, there is a unique solution p with infinite 
support. If 1 < 7 < 6/5, there are no stationary solutions with finite total mass. These 
equilibria for 6/5 < 7 < 2 are often referred to the Lane-Emden stars. 

The stability problem of the Lane-Emden star configurations has been of great interest 
both physically and mathematically. Its physical literature can be tracked back to Chan- 
dresekhar [2] and it has been conjectured by astrophysicists that the steady states for 7 < 4/3 
are unstable. To see why this conjecture makes sense, we recall the total energy E of Euler- 
Poisson system (jl.ip : 

E(p, u)= I \p\u\ 2 + -^—dx - -!- / \V<$>\ 2 dx. (1.5) 

J 2 7 - 1 07T J 

Note that E is not positive definite, which is a major difficulty of the stability question. The 
stability is related to the competition between the kinetic energy which tries to pull the gas 
away from the center and the gravitational potential energy which tries to push the gas to 
the center. The corresponding energy for the Lane-Emden equilibria of (|1.4p can be written 
in terms of the pressure integral: 

4-37 f 
E(p) = — - J pdx. 

For the explicit computation, see [I]. We see that the energy is negative for 7 > 4/3 and 
positive for 7 < 4/3 and hence one may predict the stability based on the principle of 
minimum energy. In fact, one can check that when 7 < 4/3, by constructing a scaling 
invariant family of steady states, the corresponding steady states are not minimizers of the 
energy functional (11.5j) and this indicates a possibility of certain kind of instability. 

The mathematical development of the stability theory of the Lane-Emden stars is rather 
recent. The linear stability was studied by Lin [12]: any stationary solution is stable when 
7 > 4/3 and unstable when 7 < 4/3. In accordance with the linear stability, a nonlinear 
stability for 7 > 4/3 was established by Rein in [21] using the variational approach based on 
the fact that the steady states are minimizers of an energy functional. For 7 = 4/3, in [JJ, 
Deng, Liu, Yang, and Yao showed that the energy of a steady state is zero and any small 
perturbation can make the energy positive and cause part of the system go off to infinity. 
In [5], the author proved the fully nonlinear, dynamical instability of the steady profile for 
7 = 6/5 based on the bootstrap argument and nonlinear weighted energy estimates. For 
7 = 6/5, the density of the Lane-Emden star decays rapidly without vanishing in any finite 
radius and the analysis in [5 J is not directly applied to other cases of 6/5 < 7 < 4/3 which 
attain more physical, interesting features. This is because the gas sphere for 6/5 < 7 < 4/3 
contacts with vacuum continuously, namely the density of gas vanishes at a finite radius. The 
boundary behavior of compactly supported Lane-Emden solutions is characterized as follows 
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Po( r ) ~ {R — r ) 7-1 f° r r ~ R 



(1.6) 



near the boundary. The presence of this vacuum boundary makes the problem challenging 
and very interesting and the nonlinear stability question has been an open problem for a long 
time. The purpose of this article is to establish a nonlinear instability theory of compactly 
supported Lane-Emden steady stars for 6/5 < 7 < 4/3 in the presence of the vacuum 
boundary. 

In the next section, we introduce Lagrangian coordinates to deal with the vacuum states 
satisfying (jl.6p . derive the Lane-Emden equation in Lagrangian formulation, and state the 
main result of this article. 




2 Lagrangian formulation and Main result 

2.1 Physical vacuum 

The boundary behavior given in (jl.6p naturally arises from the physical system such as the 
Euler-Poisson system or Euler equation with damping and it is referred to the so-called 
physical vacuum in the context of compressible fluids and gases [13] • In general, when the 
initial density function contains a vacuum, the vacuum boundary T is defined as 

T = cl{(t, x) : p(t, x) > 0} n cl{(t, x) : p(t, x) = 0}. 

By introducing the sound speed c 

Kjp 2 for polytropic gases), 
we call a vacuum boundary T physical if 

dc 2 

-oo< — <0 (2.1) 
on 

in a small neighborhood of the boundary, where n is the outward unit normal to T. Despite 
its physical importance, even the local existence theory of classical solutions of compressible 
Euler equations featuring the physical vacuum boundary was established only recently. This 
is because if the physical vacuum boundary condition (|2.ip is assumed, the classical theory 
of hyperbolic systems can not be directly applied. Indeed, Nishida in [20J suggested to 
consider a free boundary problem which includes this kind of singularity caused by vacuum, 
not shock wave singularity: under the physical vacuum condition (|2.ip . the boundary is 
supposed to move with a finite normal acceleration, and it becomes part of the unknown. 
Hence, one needs to deal with a moving boundary. Very recently, the local existence theory 
of the physical vacuum states for three-dimensional compressible Euler equations has been 
given independently by the author and Masmoudi [8] and by Coutand and Shkoller [3] in 
Lagrangian coordinates with different approaches. For more detail on physical vacuum, we 
refer to [9] and references therein. 

There are some existence theories available of Euler or Euler-Poisson equations in other 
smooth vacuum states when c is continuously differentiable; for instance, see |13l 1171 [TBI 119] . 
For more discussion on the existence theory concerning vacuum states, we refer to [H |9]. 
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In accordance with the recent advancement of physical vacuum [31 [8] , we will study the 
Euler-Poisson system (|l.ip or (|1.3j) in a suitable Lagrangian coordinates and formulate the 
problem in such Lagrangian coordinates. 

2.2 Derivation of the system in Lagrangian coordinates 

First we will write the Euler-Poisson system (|l.ip in Lagrangian coordinates. Let rj(t, x) be 
the position of the gas particle x at time t so that 

7] t = u(t, rj(t, x)) for t > and r/(0, x) = t]q in Q . (2.2) 

Here t]q is not necessarily the identity map but diffeomorphic to the identity map. Our choice 
of r/Q will depend on the initial density profile and in fact, in our setup, the identity map will 
correspond to the equilibrium state. 

The followings are the Lagrangian quantities: 

v(t, x) = u(t,r](t,x)) (Lagrangian velocity) 
f(t,x) = p(t,rj(t,x)) (Lagrangian density) 
i&(t, x) = <S>(t,r)(t,x)) (Lagrangian Potential) 

A = [Dr]]^ 1 (inverse of deformation tensor) 

J = det Drj ( Jacobian determinant) 

a = J A (transpose of cofactor matrix) 

We use Einstein's summation convention and the notation F,k to denote the A; th -partial 
derivative of F. In this subsection, we use i,j,k,l,r,s to denote 1, 2, 3. The Euler-Poisson 
equations (|l.ip read as follows: 

/ t + 74^ = 

fvl + KAfr, k = -fAf^, k (2.3) 

Since 

J t = JA\v\j , 
together with the equation for /, we find that 

fJ = p(0)J(0) = Pin detD Vo 

where pi n is a given initial density function. We now choose r/o so that 

fj = p in detD7] = p (2.4) 

where p is the equilibrium density profile of the Lane-Emden star given by (|1.4j) which is our 
object for the instability theory. See also Section [2.31 for the Lane-Emden equation. 

By using the relation A\ = J _1 a^, we see that the system (12. 3D is reduced to the following: 

pvj + KaHptJ^)* = -pA*% k 
A k i {A l ^, l ), k =A-KpJ- 1 



4 



along with 

4=v\ (2.6) 

Now let w be 

w = Kf- x . (2.7) 

Note that t^jW is the equilibrium enthalpy. We prefer to work with the enthalpy w to the 
density p, since w behaves like a distance function near the boundary irregardless of the 
values of 7 under the physical vacuum condition (jl.6p . See Section [231 for the basic property 
of w. We will treat w as the weight function. This viewpoint is important not only for the 
existence theory but also for our instability analysis. 

By using the Piola identity: a^, fc = 0, we see that the equation (|2.5p takes the following 
form 

w a w\ + (w 1+a Ap- l / a ), k = -w a A^, k 

A ! i(A l i ^, l ),k = ^K- a w a J- 1 

where 

a = - J- . (2.9) 
7-1 

Here a has been introduced for the notational convenience. We will use both a and 7, which 
are related through (|2.9p . in the equations and the estimates throughout the article. For 
instance, the range of the adiabatic exponents of our interest reads in terms of a as follows: 

6/5<7<4/3 3<a<5 

For our purpose of the instability theory, it is sufficient to consider the spherically sym- 
metric motion. What follows next is the derivation of the spherically symmetric Euler-Poisson 
flows in Lagrangian coordinates. 

Let 

T](t,x)=Z( k t,r)x (2.10) 
where r = \x\. Since rj t = = v, v(t,r) = r£ t . Since d{ = ^rd r , note that 

Thus J and \Drf\~ 1 can be written as 

J = t 2 (t + irr) ; [Dt,]- 1 = f~ ^ l\ rr)r (xV) (2.11) 



and hence A\ is given by 

A k = 



A\ = %- f xkx \ (2.12) 



and the gradient Afd k is given by 
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and the Laplacian A\dk{A\d{) is given by 



A*d k (A l A) = A k t d k ( d r ) = A^l 1 d r + x'Atd* ( } dr 

3 i, r r \ d r r ( d r 

I , 7 Of 



£ £(£ + &-r)/ r(£ + £ r r) £ + &r V^ + ^r 
-d r \ -r—z — I + : 7Z 2 : — r<9 r 



1 A (M-A 



(f + { r r)(£r)2 Vf + &r 
Thus the Poisson equation in (|2.8|) for spherically symmetric flows takes the following form: 

Based on ()2,10p , (|2.11[) , and (I2.12D , we see that the momentum equation in (12 .8p for spherically 
symmetric flows can be written as an equation for £: 

w^u + ^d r (w 1+a lf(Z + trr)]- 1 ^) + ™" <S r = (2.14) 
r V / r(£ + £ r r) 

for i > and < r < R where R is the radius of the Lane-Emden star. We remark that 
no boundary conditions are necessary to construct smooth solutions for (|2,14p due to the 
degenerate weights OH]. Note that from (12.13P the potential term can be also written as 



-*r = / -—W a S 2 ds 



pdx^j . 



This potential term is essentially lower-order with respect to differential structure at least for 
a short time and moreover, it has the right weight w a . Hence, together with the regularity 
of the weight w in Lemma 13.31 the local existence of smooth solutions to (|2.14j) can be easily 
established based on the theories developed by the author and Masmoudi [8] or by Coutand 
and Shkoller [3]. In this article, we take the existence of smooth solutions for granted and 
focus on the instability analysis. 

2.3 Lane-Emden star configuration in Lagrangian formulation 

In this subsection, we will identify the Lane-Emden stars satisfying (jl.4p in our Langrangian 
formulation. The static equilibria of the Euler-Poisson system under spherical symmetry 
governed by (|2.14p can be found by setting £ = 1, since pi n = p and thus J = 1. Letting w 
be the corresponding enthalpy profile (|2.7p . it is clear that w satisfies the following ordinary 
differential equation 

1 w a f r Air 

-d r (w 1+a ) + - T / — w a s 2 ds = (2.15) 

or equivalently 

2 47T 

Wrr + ~Wr + . = (2.16) 

r (1 + a)K a 
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The equation (|2.16p is the so-called Lane-Emden equation. As mentioned in the introduction, 
this equation has been studied extensively. In particular, we recall the well-known existence 
result from [21 [12]: Supplemented with the following normalized boundary conditions 

w{0) = 1 and w r {0) = 

for a given finite total mass M, there exist a ball-type solution w to the Lane-Emden equation 
(|2.16p and a finite radius R when 1 < a < 5, equivalently 6/5 < 7 < 2 such that 

(i) w > for < r < R and w(R) = 0; 

(ii) -00 < w r < for < r < R; 

(hi) w satisfies the physical vacuum condition (II. 6p or (|2.ip . 

It turns out that the Lane-Emden configuration w enjoys better regularity. In Lemma [3 .3\ 
we will study the higher regularity of w, which will be helpful to infer the suitable regularity 
of the growing mode solution. 

2.4 Main Result 

To state the main theorem, we write the equation (|2.14p in a perturbation form around the 
equilibrium state given by £ = 1 and ^ = 0. Letting £ = 1 + £ where \C\ <C 1, we obtain the 
equation for ( as follows: 



C) 2 „ , ,2,, : . , . j 

(2.17) 



W a ( tt + ^—^8 r ( W i+Q [(l + 0^1 + C + Crr)]-- 



^w a s z ds = 0. 



(1 + C) 2 r 3 7 * 

Our main result establishes the fully nonlinear dynamical instability of the Lane-Emden 
equilibrium stars satisfying (|2.16p under the Euler-Poisson dynamics (|2,17p . In the statement 
of the following theorem, for any given 5 > and 29 > 5, we define 

V 1 ln^ (2.18) 



Mo ^ 



where ^/JIo is the sharp linear growth rate. 

Theorem 2.1. Assume that 6/5 < 7 < 4/3 (3 < a < 5). Then there exist functions spaces 
Zq and Z" as well as constants 9 > and C > such that for any sufficiently small 5 > 0, 
there exists a family of solutions (dt( s , ( s ) to the Euler-Poisson system (|2.17p for t £ [0, T) 
with T > T s so that 



(d t ( s (0),( s (0)) <C5 but sup (<%C 5 (t), C 5 W) 



> 



The precise statement of Theorem 12.11 is given in Theorem 16.21 and the spaces Zq and Zf 
will be clarified in Section [H 

Remark 2.2. The escape time T s is determined through (|2.18p by the linear growth rate 
yfjM). We note that the instability occurs before the possible breakdown or any collapse of 
strong solutions. We also remark that the instability occurs in the Zq norm, which is indeed 
characterized by the instantaneous physical energy (|4.3p . 
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Remark 2.3. The above result shows that no matter how small amplitude of initial perturbed 
data is taken to be, we can find a solution such that the corresponding energy escapes at a 
time T s : there is no stabilization of the system. We conclude from this that all Lane-Emden 
steady star configurations for 6/5 < 7 < 4/3 are nonlinearly unstable. 

Remark 2.4. The instability given in Theorem \2.1\ is not related to the vorticity nor to the 
geometry of the boundary, nor to the topological change of the domain. In fact, our result 
indicates that for gases of heavier molecules 6/5 < 7 < 4/3, the self- gravitation may be strong 
enough to attract matter to the center. An interesting problem would be to investigate the 
possible gravitational collapse phenomenon rigorously. 

Up to our knowledge, this work is the first rigorous result to address the nonlinear in- 
stability of the Lane-Emden stars for 6/5 < 7 < 4/3 under the Euler-Poisson system. Our 
results together with the earlier works [U El [21] now complete a nonlinear stability and in- 
stability theory of the Euler-Poisson system for the Lane-Emden equilibrium stars, which 
provides a rigorous account of the longstanding conjecture by astrophysicists [2]. Also, we 
believe that the above result for the Lane-Emden stars and the methods developed in this 
paper will shed some light on other stability problems in the context of free boundary. 

We mention that based on strong solutions in the mass Lagrangian coordinates [6], the 
nonlinear instability theory for the Lane-Emden stars under the Navier-Stokes-Poisson system 
has been recently established by the author and Tice in [TO]: it was shown that there's 
no stabilization of the system for 6/5 < 7 < 4/3 in the presence of viscosity around the 
Lane-Emden equilibria. On the other hand, away from the Lane-Emden equilibria, there are 
interesting works on the stability theory of the Euler-Poisson system (jl.ip by Luo and Smoller 
[15} I16j. where a nonlinear dynamical stability of rotating stars for 7 > 4/3 was established 
based on a variational approach by the energy-Casimir technique and a finite time stability 
of the entropy weak solutions was proven. 

The proof of Theorem 12. II is based on a bootstrap argument from a linear instability to a 
nonlinear dynamical instability. Passing from a linearized instability to nonlinear instability 
demands much effort in the PDE context since the spectrum of the linear part is fairly 
complicated and the unboundedness of the nonlinear part usually yields a loss in derivatives. 
The key is to identify the right function spaces where the nonlinear instability would occur 
and this requires a careful nonlinear analysis. We develop nonlinear energy estimates for the 
whole system so that: first, the nonlinear estimates can be closed; and second, their interplay 
with the linear analysis can complete the argument. For this particular problem, due to the 
degeneracy caused by the physical vacuum, the novel nonlinear weighted energy estimates 
together with Hardy inequalities are employed. 

A key step is to derive Proposition 16. 1\ which is an important estimate for the bootstrap 
argument. There are a few new ingredients. In order to avoid the coordinate singularity 
at the origin due to the consideration of spherically symmetric flows as well as the vacuum 
singularity due to the physical vacuum boundary, we will apply the time differentiations 
to obtain the nonlinear energy inequalities systematically and then recover all the spatial 
derivatives from the equation by exploiting the elliptic estimates. A pivotal tool to justify 
such a process is the Hardy inequality; in particular, we will make use of the localized Hardy 
inequalities near the origin and near the boundary and such Hardy inequalities allow us to 
overcome the coordinate singularity near the origin as well as the vacuum singularity. In 
addition to the homogeneous weighted energy £ , we further introduce the nonlinear energy 
(£. This is because £ works efficiently for the bootstrap argument from the linear analysis to 
the nonlinear solutions, while the nonlinear energy (5; captures better the intriguing structure 
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of nonlinear, degenerate pressure gradient term. We will derive the energy inequality for (£ 
having the form of 

— < B£ + lower-order derivative terms 
at 

where f3 is no greater than the sharp linear growth rate ^/Jiq. To do so, we identify the 
quantities whose L°° bounds can be controlled within the total energy and use them crucially 
in controlling the nonlinear terms. With the equivalence of two energies E and (£, we can 
close the argument. 

The rest of the paper proceeds as follows. In Section [3l we construct the largest growing 
mode solution to the linearized Euler-Poisson system, study the regularity of such growing 
mode profile, and show that our growing mode solution grows at the largest possible rate. In 
Section HI the instant temporal energy and the total energy are introduced and we establish 
the equivalence of two energy norms via various Hardy inequalities. In Section [5l we derive 
high-order nonlinear weighted energy inequalities for the nonlinear energy norms and show 
that the nonlinear energy is equivalent to the homogeneous total energy. Based on the linear 
growth and the nonlinear estimates, we then prove the bootstrap argument and Theorem 12. II 
in Section El 

Throughout the paper, we will write C for a generic constant unless the distinction be- 
tween constants is vital. And we use 21 ^ *B to denote that 21 < C*B for a constant C > 0. 



3 The largest growing mode and Linear instability 

In this section, we study the linearized Euler-Poisson system around the Lane-Emden stars 
and show the linear instability of the Lane-Emden stars whose indices are 6/5 < 7 < 4/3 
(3 < a < 5) in the sense of Lin's criterion [12J. We will first derive the linearized equation of 
(|2.17p . Notice that by Taylor's theorem, for sufficiently small £, the nonlinear pressure term 
in (|2.17p can be written as 

_ l+q 1 _|_ Q/ 

[(1 + C) 2 (l + C + Crr)] * =1 (3C + Crr) + /»(C,Crr) (3.1) 

where h is a smooth function in both arguments £ and £ r r and h(C,Cr r ) = 0(\(\ 2 + |£ r r| 2 ), 
and thus the linearized equation of (|2.17p reads as 

w«Cu + l^llC " — [w 1+a (SC + 0)] r = (3.2) 
r ar r 

where we have used (|2.15p . 

To find a growing mode, we seek a solution to (|3.2p of the form r) = e xt (f>(r). Then cf> 
would satisfy 

\ 2 w a 4> + 4[W +a] ' cP - i±^[ w 1+a (30 + ct>'r)}' = (3.3) 
r ar 

V v ' 

(*) 

where we denote ' = The (*) in (|3.3p can be rewritten as 

(*) = (4 - 3^)1^0 - i±^V+>V 4 ]'. 

a r ar q 
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Multiply (|33J) by r 4 to get 



AW 4 = ^±[ w l+<*rU'}' - (4 - 3^^)[w 1+a }' r 3 </> (3.4) 
a a 

We denote the right-hand side of (|3.4h by L<p: 

Lcf) = [w 1+l Vy ]' - (4 - 3^±^) r 3 (3.5) 

Then L is self- adjoint and hence A 2 is real. 

We recall that in [12], the stability criterion was introduced based on the eigenvalues: 
w a (~ p) is called neutrally stable if A 2 < for all eigenvalues A and unstable if Af > for 
some eigenvalue Ai; and that it was shown that w a (~ p) is unstable for any 3 < a < 5 
(6/5 < 7 < 4/3) in mass Lagrangian framework. In the next subsection, we verify Lin's 
linear instability result in our Lagrangian framework by showing the existence of the positive 
largest eigenvalue for 3 < a < 5. 



3.1 The existence of the largest growing mode 

We start from Lcp = pw a r 4 cp where L is defined in (13. 5D and \i = A 2 . It is well-known that 
the largest eigenvalue po of this equation can be found by the variational formula 

no = su p|j^ : Q(4>) < oo, 1(0) < ooj 



where 



Q(cf>) = (Lcf>,cf>) = -^-^- f R w 1+a r\(j)') 2 dr-(4-3^-^) f [w 1+a ]' r 3 cf) 2 dr, 
a Jo a Jo 



I(cf>) = (w a r A cf > , cj>) = I w a r 4 cf> 2 dr. 



o 



Here we have used ( , ) to denote the standard one-dimensional inner product on (0, R). 
Define a norm for any cp E C°°([0,R]), 



" Ft pR 
,1+Q!„4/ ±l\2 j„ / r„„l+an' 3 ,2 i , / n „a4,2, 



w L+a r\cj)'ydr- / [ti/ +a ]'r'V / dr + / w a r*fdr. 
Jo Jo 

Let -ff = C°°([0, i?]) in the above norm. Since jjjT = f° r an y nonzero constant c, the 

variational problem can be restated as finding the maximum [iq of Q(4>) on H under the 
normalization I(cf>) = 1. The next proposition shows that the supremum jiq can be achieved 
on H. 

Proposition 3.1. There exists a (fio 6 H such that I((/>o) = 1 and Q(4>o) = Ho > 0. 

Proof. The first part in Q is negative and the second part in Q is positive because w' < 
and 4 — 3^j— > since a > 3. Thus we first want to make sure that hq > 0. This can be 
easily checked by taking = 1: 

Q(l) _ Q - (4 - 3^) /* K+«]'rW 
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Note that the positive part (the second term) of Q is uniformly bounded by I because from 
(I2,15p . we have that — ^^[w 1+a ]' = \ f$ j^w a s 2 ds is bounded. This implies that /xo is 
finite. To verify that //o is attained on H, we will need a compactness argument. Let (f) n be 
a maximizing sequence so that 

Q{4>n) f~ Mo as n — * 00 an d I{4> n ) = 1 for all n. 
Let i^o be its weak limit. By the lower semicontinuity of weak convergence, we have 

r-R r-R 

liminf / w l+a r 4 (<j)' n ) 2 dr > / w 1+a r A (<p' ) 2 dr 
Jo Jo 

R rR 



liminf / w a r (4> n ) dr > w a r (c/>o) dr 
Jo Jo 

The next step is to prove the following compactness result. 

Claim 3.2. The mass is not lost in the limit: I(4>o) = 1. Moreover, we have the compactness 
result for the positive part of Q: there exists a subsequence {4>n k } of {4> n } such that 

r-R r-R 

/ -[w 1+a ]' r 3 <f>l k dr ->■ / -[w 1+a ]' r 3 <f) 2 dr as 00. (3.6) 
Jo Jo 

Since the above claim implies Q((fto) = Ho and / (^0) = 1> the conclusion of the proposition 
follows from the claim. 

The first claim can be shown by a standard scaling argument. Suppose I((fio) = /3 2 < 1. 
Then we would have 

J (-g) = 1 and Qi-J) = lp > 

which contradicts to the definition of /iq. Hence, I((j)o) = 1. Next, in order to see (j3.6H . we 
rewrite it by using (|2.15[) 

* - [w i+"}> r^ldr = £ (^£^- W a s 2 ds) vf^fidr 

Then since -3 J* Q r jJ^w a s 2 ds is uniformly bounded from below and above on (0, R), the positive 
part of Q is bounded from below and above by /: 

R rR / 1 rr 4- \ r R 

-[w 1+a ]'r i <p 2 n dr = J U J — w a s 2 ds)w a r 4 <t> 2 n dr~ ^ w a r^ 2 n dr = I(J> n ) 



Now since I{4> n ) — > !{4>o) with both being one, (|3.6p follows. □ 

It is straightforward to check that ^0 is in fact an eigenfunction corresponding to fio for 
(|3.5|) : LipQ = ^ow a r 4 (j)o and therefore, C(^ r ) = ev / w*0 o (r) is a largest growing mode solution 
to the linearized Euler-Poisson equation (|3.2p . 
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3.2 The behavior of </>o near the origin and near the boundary 

We record the second-order ordinary differential equation which the growing mode 4>$ satisfies: 

» w a rUo = ^[w 1+a r%]' - (4 - 3^)[w 1+a ]> r 3 </> (3.7) 
a a 

We can further deduce the regularity of 4>o from (|3.7p based on the behavior of the Lane- 
Emden solution w. To do so, we first derive the higher regularity and integrability of w. 

Lemma 3.3 (Regularity of w). Let 1 < a < 5 be given. Let w be a ball-type solution to the 
Lane-Emden equation (|2.16p . Then 

1. w is analytic near the origin. Moreover, 

w (r) = 1 - br 2 + 0(r 4 ), r~0 

for some positive constant b > 0. Also, (d 2k+1 w)(0) = for any non-negative integer 
k > 0. 

k 1 

2. d l r w is uniformly bounded on (0,R) for each < i < a + 2 and also u> _ 2 _ d k+1 w is 
uniformly bounded on (0, R) for each 1 < k < 2a + 1. /n addition, w enjoys the 
following integral regularity. 

R 

w a+j r A \di +l w\ 2 dr < oo 



for each < j < 3a + 3. 

The proof is elementary, but it has not been found in the literature, so we will provide a 
brief argument. 

Proof. We start with the first item. Write the Lane-Emden equation (I2.16P as 

rw rr + 2uv + crw a = (3.8) 

where c = ri+a)K a > i s a constant. After dividing (|3.8|) by r and taking the limit r — > + , 
we immediately deduce that 3u; rr (0) + c = 0. Take d r of fj3.8j) : 

™ rrr + 3w rr + c(arw a_1 w r + w a ) = (3.9) 

By dividing it by r and using the fact that 3w; rr (0) + c = 0, 

3w rr - 3w rr (0) , „_! ^"-^"(O) 
w rrr . H h c(aio w r H ) = 

and by passing to the limit r — > + , we see that 4w rrr (0) = 0. Take one more derivative 

rwrrrr + 4u; r rr + c(axu> Q-1 uv r + 2aw a_1 w r +a(a — l)rw a_2 (tf r ) 2 ) = (3.10) 

=— acrw 2a ~ 1 

and we deduce that 5w rrrr (0) — ac 2 = 0. More generally, we obtain the following. 
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Claim 3.4. For any k > 1, the solutions to (|3.8p satisfy the following 

rd 2k w + 2kdf- x w + td 2 r k - 2 {rw a ) = (3.11) 

for < r < R such that 

(a) d 2k - 2 (rw a ){0) = and d 2k (rw a ){0) = 

(b) (d 2k ~ l w){o) = o. 

We will prove the claim by induction on k. First note that (a) and (6) are satisfied for 
k = 1 and A; = 2 by the above computations. Now let us assume that it is true for k = n. 
Then by dividing (|3.1ip when k = n by r and taking the limit r — > + , we deduce that 

(2n + l)(d 2n u;)(0) + Z n = (3.12) 

where l n = lim r .^o+ [c<9 2n_2 (™ a )/r]. Now we formally take d r of (|3. 1 1 () : 

rd 2n+1 w + (2n + l)d 2n w + cd 2 "" 1 ^™") = 

By using (|3.12j) . we see that 

^ + ( 2 „ + 1) 9>-(9»(0) + <-K)-'. _ 
Take the limit to get 

(2n + 2)(5 2n+1 u;)(0) + ca r 2n (ru; a )(0) = 

Since 3 2n (r'u; Q )(0) = 0, we deduce that (d 2n+1 w)(0) = 0. Also, the direct computation shows 
that 9 2n+2 (ru> Q )(0) = 0. And thus (a) and (b) hold for k = n + 1. Hence, we deduce that 
it; is infinitely differentiable at r = and the above estimation on the derivatives implies the 
analyticity of w near the origin. This finishes the above claim as well as the first part of the 
lemma. 

In order to prove the point 2 of the lemma, it suffices to look at w near r = R since 
w is regular near the origin. Recall that w r = 0(1) for r ~ R. From (j3.8|) . we see that 
w rr = 0(1) + 0(w a ) for r ~ R. Similarly, from (HUD and (|3TT0|K w rrr = 0(1) + O(u)^ 1 ) and 
w rrrr = O(l) + 0{w a ~ 2 ). Inductively, we obtain \&lw\ = 0(1) + 0(w a ~ j+2 ). Note that dlw 
is not bounded any longer if j > a + 2. However, the integral 



/ w a+j r A \di +1 w\ 2 dr <C + C 
Jo Jo 



w 



stays bounded as long as j < 3a + 3. On the other hand, if we consider w~^~ d k+1 w, 
then we see that \w^~ d k+1 w\ = 0(1) + 0(w a ~ k+1J[ ~), which is uniformly bounded if 
a — k + 1 + >0=^1<A;<2q + 1. This completes the proof of the lemma. □ 

The following lemma concerns the behavior of (f>Q near the origin. 

Lemma 3.5. Let (fro £ H be the solution to ()3.7|) . (j>o is analytic at r = and moreover, 
= a + 0(r) around the origin where a is a constant. 



13 



Proof. We rewrite (|3.7p in the following form 



+ 



4 (1 + a)w' 
r w 



a(4-3 ) — + — 

a rw 1 + aw 







Let P(r) and Q(r) be coefficients of (fr' and (fro respectively. Then by Lemma 13.31 i"P(r) and 
r 2 Q(r) are analytic at r = 0, which implies r = is a regular singular point. Therefore, we 
may employ the Frobenius method [1] (p. 215). Let po and go be the zeroth order term of 
rP(r) and r 2 Q(r) respectively. It is easy to see that po = 4 and go = 0. The indicial equation 
r(r — 1) +Po r + Qo = has two roots r\ = and r2 = —3. Hence, by the Frobenius theorem, (fro 
has a power series representation of y\{r) = Yln=o a « r?1 or 2/2 ( 
However, yi (r) is impossible since (fro £ H and hence </ 



oo 



y x (r) In r + r 3 £ n=0 Un 



o„r" 



□ 



Next, we show that 0o enjoys the better integrability near the boundary than the inte- 
grability dictated by H in Proposition 13.11 This will be done through the Hardy inequality 
and a bootstrapping argument through the equation. In addition, we show that <pQ indeed 
belongs to some weighted Sobolev spaces. 



Lemma 3.6. Let 4>o be the solution to (I3.7p . Then 

1. (fro has the following integrability: for any < f3 < a, 



[ w a - p r i (frldr+ [ 
Jo Jo 

Moreover, for any z > 1, 



w 1+a -Pr 4 (4>' ) 2 dr <oo. 



w 



z-^cfrl 



dr < oo. 



o 



2. (fro has the following regularity: for 1 < k < 2a + 1, 



w 1+a+k r 4 (d^ +1 (fro) 2 dr < 



oo. 







Proof. We start with the first item. j3 = corresponds to the regularity of H in which (fro 
was constructed. From the regularity of H , and by using the Hardy inequality - see Lemma 
14.21 as well as (|4.1ip - we deduce that 



R 



w a ~\^ 2 



\dv< [ R w a+l r 4 (fr 2 odr + T 
Jo Jo 



w a+1 r\(fr'ofdr < oo. 



In order to bootstrap this better integrability for (fro into the one for (fr' , we will make use 
of the elliptic structure of the equation (|3.7j) . Multiply (|3.7j) by w^ 1 ^ and integrate it over 
(0,1?): 



1 + a. 



\w 



a 



w 



dr + 



R 



HoW a 1 r A (fr 2 ) dr 



(*) 



R 



-(4-3 



1 + a [w 1+a ]' 



(3.13) 



r^cfr^dr 



a 



IV 
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Note that by Lemma [3] 



(**) 



R 



w a 1 r 4 ip^dr < oo. 



For the (*), by integration by parts, we obtain 



a 



1 + a 



From (|2.16p . we see that 



R 



(*) = / w a r\<t>' Q ) 2 dr - / w'w^r^atfodr 



R 



w a r 4 ((P' Q ) 2 dr + 



a — 1 



+ 



w"w a - 1 r 4 ((l)o 



(w') 2 w a - 2 rU 2 dr 

r R 

) 2 dr + 2 / v/vp-^ffidr 
Jo 



(*)i 



(3.14) 



l(*) 



R 



w^rUldr 



R 



rv—1 4 il i 

w r <p dr < oo 



and hence, together with (|3.14p from (|3,13p we deduce that 



R 



w a r A {4 ) ' Q ) 2 dr 



a — I 



R 



{w') 2 w a - 2 r 4 (f)ldr < oo. 



This establishes the case of f3 = 1. The general case can be done by a similar bootstrapping 
argument: assume that w^^r^c^dr < oo for < f3 < a. By multiplying (|3.7p by w~&'- 
and integrating to derive that 



R 



w a+1 -Pr\0' o ) 2 dr + 



a-0 



(w') 2 w a -P- l rU 2 dr 



is bounded. We remark that the condition on f3 is necessary to keep the second integral stay 
non-negative; otherwise, we won't be to close the argument. Now by the Hardy inequality 
again in Lemma 14.21 we conclude that w z ~ 2 r 4: (j) 2 ) dr is bounded as long as z > 1. 

We move onto the item 2 regarding the higher regularity of 0o- The key to achieve this 
higher regularity for </>o is to capture the right behavior of L with respect to the spatial 
derivatives. We will follow the idea in [8]. We first write (|3.7|) as follows: 



wcj)Q + (1 + a Wo + — 4>' 



a(4 - 3— — )— 0o + — — < 
1 + a r 1 + a 



(3.15) 



Take the derivative of (|3. 15[) and write it as 



+ (2 + Q»o' + ^<t>0 



(l + aKV + (^)Vo 

a w' a^o 
a(4 - 3 )— (po + 



1 + a r 



l + a 



(3.16) 
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so that the left-hand-side of (|3. 16[) keeps the self-adjoint structure with different weights: 



Aw 



w(f% + (2 + o>Vo + —<t>o 



Hence, by multiplying (|3.16p by w 1+a r 4 4>' and integrating, we get 



R 



f R dm 

,»fdr= / [(l + a )w" + (— )> 1+t V 4 (<#,) 2 dr 
J o r 



R 



«( 4 - 3— — — 0o + 7- — < 
1 + a r 1 + a 



w 1+a r 4 (j) dr 



It is easy to see that the right-hand-side is bounded because of the property of w and since 
00 £ H. The higher-order estimates can be inductively done in the same manner. By taking 
more derivatives of (13. 16ft . one can derive that 



Aw 



w 



d r fc+2 o + (k + l + a)w'd* +1 cp + — d k r +l 



T 



(3.17) 



where T consists of lower-order terms such as (d k+l w)4>' , (d k w)<pQ, w"d k 4>o and so on. The 
left-hand-side of (13. 17ft is written as 



Aw 



wd k+2 <p + (k + 1 + a)w'd k+1 (/) + —d k+l 



w k+a r 4 



W 



k+l+a r 4 d k+l ( 



and thus by multiplying by w +a r dScfio and integrating, we get 



w 1+a+k r A \d k+l 4> \ 2 dr = / T -w k+a r A d k (t> Q dr 



What remains to show is that the right-hand-side is bounded. Suppose that T has the form 
of {di +1 w)(d k+1 ~ " J (fto) for 1 < j < k; other terms are even lower than these terms and they 
can be treated in the same way. Then 



R 



T ■ w k+a r 4 d k <f) dr => / (w^rd} +l w)(w 



R 



1±1 ^r 2 d k+1 -^o)(w^r 2 d k ( Po)dr 



and thus by using Lemma 13.31 and by using Cauchy-Swartz inequality, we deduce that it is 
bounded. Hence, we conclude that w 1+a+k r 4 \d k+1 <f)Q\ 2 dr < oo as long as 1 < k < 2a + 1. 



This completes the proof the lemma. 



□ 



Remark 3.7. Lemma \3. 5\ and \3. 6\ show that the growing mode 4>q does not display an erratic 
nor singular behavior near the origin and the boundary. Indeed, this growing mode profile 
belongs to the function spaces of our interest, namely has a finite total initial energy for 
3 < a < 5; see (|4.5p and (|4.6p . And hence, we can use it as an initial perturbation for the 
nonlinear problem. 
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3.3 The linear growth rate 

We will now show that y/JiQ is the dominating growth rate for the linearized Euler-Poisson 
equation. We begin by writing the linearized equation (|3.2p for = £ 



w a r 4 ^ tt = (3.18) 

where L is the self-adjoint linear operator given in (13. 5p . Here, instead of £, we use ^ to 
distinguish the linear analysis from the forthcoming nonlinear analysis and we note that 
= ^(t,r) while the <j> in (I3.4p and (I3.5D depends only on r. 
We now introduce the following weighted norms: 



2 X'= [ R w a r 4 fdr, 
Jo 

2 + « f R „„a+l 4 ,2, 



= 

a Jo 



(3.19) 



w a+1 r*tfdr. 



\x a - 



The first result is on the estimate of \& satisfying f|3. 18[) in 
Lemma 3.8. For every solution \l/ to f)3. 18|) i/iere exisi Cu ,Cn j > suc/i i/iai 

||*(t)||x-, ||*t(i)||x« < C Mo e^(||*i(0)||x« + ||*(0)||x- + ll*(0)|M (3.20) 
and /or any i > 1, 

< C Mii e^(||vI/ t (0) + ||*(0)||x« + ||*(0)[|y«) 

* (3 21) 

+C /iO)i ^(||^ +1 *(0)||^ + ||W0)||y«). 

i=i 

Proof. Take the inner product of (|3.18p with vpj. Then we get 

(w a r 4 * tt ,%) = (L*,%) <=> ^(™ a r 4 ¥ t ,¥ t ) = ^(Ltf,tt). 

The above equivalence comes from the self-adjointness of L. Next integrate the above in time 
t to get 

(w a r 4 %(t), ¥ t (t)) = (L*(t), ¥(t)) + KV% t (0), tf t (0)) - (L*(0), #(0)). (3.22) 

Since (L^ (t) , ^ (t)) < fio(w a r 4l ^(t), 1 ^(t)) for all t because of the definition of /io and since 
-(Ltf(0), tf(0)) < ||*(0)|| yce , from (j332|) we see that 

||*t(t)||^a < Mo||*(t)||^a + ||* t (0)|&„ + ||¥(0)|&„. (3.23) 

Now since 

||*(*)||x<*< f \\MT)\\x*dT+\\V(0)\\ x >> 
Jo 

plugging this into (|3.23p . we get 

||*t(t)||x« <VJ^ f \\MT)\\x*dT + C(\\%(Q)\\ X c, + ||¥(0)||x- + ||*(0)||y«). 
Jo 
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By the Gronwall inequality, we obtain 

||*t(t)||x« < Ce^ mt (\\^ t (0)\\ x « + ||*(0)||x«» + ||*(0)||y«) and in success 
\Mt)\\ X a < CeV^*(||* t (0)|| x « + ||^(0)||x« + ||*(0)||y-). 

Note that C only depends on //q. For higher derivatives, take df of f)3. 18|) : w a r^d\^tt = Ldfo. 
Take the inner product of this with dtd\^ to get 

= fljtf(i)) + Kr*ff t (0),ff t (0)) - $¥(0)) 

< toHmtc- + ll^ +1 *(0)ll^ + IIW)|| y «. 
Hence (|3T2TD follows. □ 



The estimates of ¥ in Y a space can be derived by using the self-adjointness of L again 
and by Lemma 13.81 



Lemma 3.9. For every solution ¥ to (|3.18p there exists C„ and > such that 
\Mt)\\h < C, o e 2 ^\\\M0)\\x« + [|*(0)|&. + ||¥(0)|£.) 

and for i>l, 

II W)IIy« < Cn,i^(\\%{Q)\\x>> + ll*(0)l&- + ll*(0)|&-) 

i 

+c^Yl(\\% +1 m\\x« + \\%meY>>)- 

3=1 

Proof. Here we give a proof for i = 0. Other cases of i > 1 can be treated in the same way. 
By using (]3.5p , we rewrite the equation (|3.22p as 

\Mt)\\Y« = - (4-3^) [ R [w 1+a ]'r 3 ^dr + [|¥ t (0)|&« - (L*(0), ¥(0)). 

Notice that since — — - - ^ u; a from (|2.15p . the second term in the right-hand-side is 
bounded by C \\ V(t)\\ Xa . Moreover, by recalling that -(L¥(0), ¥(0)) < ||*(0)|| y „, we deduce 
that 

[|*(*)[|^ < C\\9(t)\\x- + ll**(0)[|i- + ll*(0)||y«- 
And therefore, from Lemma 13.81 the conclusion follows. □ 



In the next two sections, we develop the nonlinear theory for solutions to the Euler-Poisson 
equation in the spherically symmetric motion (12.17p . 



4 The instant energy and total energy 

In this section, we will introduce the various energies and establish the equivalence of the 
temporal instant energy and the total energy for the solutions to the Euler-Poisson equation. 
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We first rewrite the equation (|2.17p as 

———2 w + r(r r + 3 4 )w r $ (r C ^ 

(1 + C) a a (4.1) 

+ r 3 (u; 1+Q /i(C, Crr)) r - w a r 4 $(r)/(C) = 
where $(r) is the prescribed function defined by 

*(r) = - <^ = -(1 + a A (4.2) 

and /(C) is given by 

f(n= (i + Q 4 -i-4C(i + C) 4 
AW (i + C) 4 

and thus /(C) = 0(|C| 2 ) for small C, and finally h is given in f)3.1|) . 

We introduce the following instant energies and the total energy for the solutions to the 
nonlinear Euler-Poisson equation (]4.ip . We denote the zeroth-order kinetic energy plus the 
lower-order potential energy for C by £ °: 

rR 1 i „, /--R 

„4|a|2 



f o_ /" ^ r 4|G| 2 + ^u; 1+a r 4 |C,| 2 dr+ / 
Jo a J 



w a r*\(\ z dr (4.3) 



The higher-order (temporal) instant energy is denoted by £ 3 for j > 1: 

£J = [ R w a r 4 \di( t \ 2 + i^w 1+a r 4 |^Cr ! 2 ^ 
Jo « 

and the total instant energy by £(t): 

N+3 

i=i 

where [a] is a ceiling function, namely [a] = min{n S Z : a < n}. Notice that the instant 
energies can be expressed in terms of X a and Y a that were introduced earlier in (|3.19p as 
follows: 



^IICtll^ + IICII^ + IICHx* 

£ j = Wdittfx* + ll^Cllf- 

For the mixed derivatives, we introduce the following notation: 

S j,k = [* w °+*r 4 \di- k d*C t \ 2 + ^^w l+a+k r 4 \d{~ k d^Cr\ 2 dr. 
Jo 1 1 a 



(4.5) 



We remark that this change of weights for the spatial derivatives in the energy norms is not a 
coincidence but rather a nature of the physical vacuum state of polytropic gases; this nature 
was discovered and analyzed in [8] for general a. The total energy including both temporal 
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and spatial derivatives is defined by 

\a]+3 j 

e(t) = £°(t)+ £ Y. £j,k - ( 4 - 6 ) 

j=l k=0 

It is clear that 

£{t) < £(t). 

We want to show that the converse holds for the solutions to (|4.ip in such a way: £ is bounded 
by £ under the following smallness assumption: 

[I2l±3 ]+1 [M±3] 

|C| + |Cr|+ £ \ w ^%<\ + £ h^^Cr|<0i (4.7) 

q=l q=l 

where 6\ is a sufficiently small, fixed constant. The validity of this assumption within the 
total energy £ will be justified in Lemma 14.71 

Our main goal in this section is to prove the following proposition, which concerns the 
equivalence of £{t) and £{t) under the smallness assumption (|4,7p . 

Proposition 4.1. Let (Cr,() be a solution to (14. lj) and moreover assume that (|4.7p holds for 
< t < T. Then there exists a constant C > such that for < t < T 

£(t) < C£(t). 

In order to prove Proposition 14.11 we first recall standard Hardy inequalities, and embed- 
ding inequalities, and derive more of Hardy type inequalities which can be adapted to our 
energy spaces induced by £. The proof of Proposition 14. II will be given in Section [431 

4.1 Hardy inequalities 

The following is the most well-known Hardy inequality applied to the Sobolev spaces: 

oo 

\g'{x)\ 2 dx. (4.8) 

We will also make use of the more general version of the Hardy inequality: 

Lemma 4.2 (Hardy inequality). Let k be a given real number. And let g be a function 
satisfying J s k (g 2 + g' 2 )ds < oo. 

1. Then if k > 1, then we have 



roc 




dx £ ! 


/o 


X 


J 



[ 1 s k - 2 g 2 ds<C I' s k (g 2 + \g'\ 2 )ds. 
Jo Jo 



2. And if k < 1, then g has a trace at x = and moreover, 



C s k - 2 (g-g(0)) 2 ds<C C s k \g' 
Jo Jo 



2 ds. 
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For the proof of Lemma 14.21 we refer to [llj. 

As an application of the standard Hardy inequality (14. 8p . we first obtain the following 
Hardy inequality localized near the origin. 

Lemma 4.3. Let u G X a n Y a be given where X a and Y a are defined in (|3.19p and a is any 
real number. Let c be a positive fixed number such that < c < 2c < R. Then there exists a 
constant C > independent of u such that 

rc / p2c r2c \ 

J r 2 \u\ 2 dr<c(j r 4 \u r \ 2 dr + J r 4 \u\ 2 dr) . (4.9) 

Note that the right-hand-side of (j4.9|) is bounded by ||w||^ a and ||u||y-o since the weight 
w is bounded from above and below from zero on (0,2c). 

Proof. Suppose u is a smooth function. Choose a smooth cutoff function x € C°°(0, oo) such 
that < % < 1, % = 1 on (0, c), and x = for r > 2c. Then 



I 



2c f-oo r 2 

2\ |2 j / I ' A u i2 j 

r x M dr = etr 



r 

(■CX3 



/■OO 

< C / I (r 2 x"u)r| 2 dr by Hardy inequality (14. 8p 

JO 



r>2c 

C / (r 2 (xu) r + 2rxu) 2 dr 
o 



/•2c |-2c /-2c 

C( / r 4 | (xn) r | 2 dr + 4 / r 3 xu(xu) r dr + 4 / r 2 (xn) 2 dr) 
io Jo Jo 

r2c r-2c 

C( / r 4 | (x^)r| 2 rf r — 2 / r 2 (x u ) 2 dr) by integration by parts 



which yields that 



2c /-2c 

r2 |x w | 2( ^ r — C / t 4 | (xt/)^ 1 2 <ir. 



Hence, we obtain 



f 

Jo 



-2c i-2c 

r 2 \u\ 2 dr < j r 2 \xu\ 2 dr < C j r 4 \(xu) r \ 2 dr 

2c r 2c 



r 4 x 2u rdf + j r^Xr^dr^j 

'■ r 2c \ 

r \u r \ dr + r \u\ drj. 



< C 

< C 

The inequality (|4.9p is easily extended to u G X a n l" a by the density argument. □ 

The term as in the left-hand-side of (|4.9p . which has the stronger weight near the origin, 
will appear in the subsequent nonlinear estimates due to the choice of the spherical coordi- 
nates in our analysis. Lemma 14.31 asserts that such a coordinate singularity can be resolved 
in our energy spaces by means of Hardy inequality. We can apply the same argument for 

-R / r-R r-R 

\x v \ 2 dr<C[ / r 2 \xv r \ 2 dr + / r z \xv\ z dr 



I 

Jo 
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and then apply (|4.9p to both terms to obtain the following bootstrapping estimates near the 
origin: 



R / fR fR pR 

2 J„ / I Jl ,„. |2j„ , / „4|, .„. |2 7 , / 4| 



dr < C r \xv r r\ dr + J r \xv r \ dr + J r \xv\ dr J . (4-10) 

Similarly, by using the general Hardy inequality in Lemma [4.2[ we can derive the following 
Hardy inequality localized near the boundary: 

R / pR i-R \ 

w a - 2 \ipv\ 2 dr <C (J w a \^v r \ 2 dr + J w a \ipv\ 2 dr\ (4.11) 

where ip is the smooth cutoff function satisfying 0<ip<l,ifj = for r < R — 2c and ip = 1 
for r > R — c. 

Next we recall the Hardy type embedding for the weighted Sobolev spaces which will be 
importantly used for the nonlinear energy estimates. Let be a smooth domain in W 1 and 
let d = d{x) = dist(x, be a distance function to the boundary. For any given a > and 
given nonnegative integer 6, we define the weighted Sobolev spaces H a ' b (£~l) by 

H a > b (n) = {d^u G L 2 (n) : [ d a \V k u\ 2 dx < oo , < k < b} 

Jn 



with the norm 

II u I 



b 

2 Ha , b = Y^ / d a \V k u\ 2 dx. 

1 — n " " 



fc=0" 

u.b 



We denote the standard Sobolev spaces by H s . Then for b > a/2, the weighted spaces H' 
satisfy the following Hardy type embedding [11] : 



H a '\fl) tf 6 -i 

with the estimate 

||u|| H 6-a/a ^ IMIff<a- (4.12) 

For instance, the inequality (|4.10p can be obtained from the localization (by the cutoff tech- 
nique as in the proof of Lemma l4.3f) of (I4.12p near the origin by taking Q = (0,R), a = 4 
and 6 = 2. 

We now turn our attention to the degeneracy near the boundary r = R. We would like 
to adapt the inequality (|4.12p to our energy spaces (|4.6p . Since there is loss of the weight w 
for each spatial derivative, the space for (Ct,C) generated by Ylk=o^ k ' k near ^ ne boundary is 
equivalent to the localization of r\^ =0 H a+k,k x fj a + k + 1 > k + l near the boundary with the above 
notation. Notice that r)% =0 H a+k+1 < k+1 C H a+n+1 < n+1 . Therefore, it is desirable to derive 
some Hardy inequalities by using these spaces: H@ +l ' 1 for a given positive number /? > 
and positive integers I > 0. First, by applying the Hardy embedding inequality ([4.12]) for 
a = (3 + I and b = I, we obtain 

II^MI »=£ 3 H^HIf^+M (4.13) 

£1 2 

where is the smooth cutoff function satisfying 0<ip<l,ip = {) for r < R — 2c and ip = 1 
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for r > R — c. Hence, by letting I = [/?] , we obtain the following: 

\M 2 dr^y2 wP +m \d*(ipu)\ 2 dr<J2 w^ +k \d k (tp 
k=o Jo k=o Jo 



u)\ 2 dr 



If > 2, we can combine with the Hardy inequality near the origin (|4.10p to derive the 
following: 

w^ +k r 4 \d k u\ 2 dr. 



I u 2 dr < CV [ 
Jo k=Q Jo 



Moreover, if we choose / = + 2m for m > in (I4.13P and if we combine with the Hardy 
inequality near the origin (|4.10p when + 2m > m + 2, we further obtain 



I H m 



r/31+2m R 

< C V / w p+k r A \d k u\ 2 dr. (4.14) 



k=0 

We have established the following. 
Lemma 4.4. Let d k u G 

x P+k f or each < k < \/3] + 2m w/iere f/3] + m > 2. TTien 
u G H m (0,R) and moreover, the inequality (I4.14p holds. Similar results are valid for Y@ : if 
d k u G Y^ + * /or eac/i < k < \p\ + 1 + 2m, ffren u G # w (0, 12). 

As a result of Lemma |4.4( we get the following L°° embedding of our energy spaces: 

Lemma 4.5. Let d k u G X /3+fc for each < A; < + 2 w/iere [/?] > 1. Then u G L°°(0, 12) 
i/te following estimate: 



l«IIOO 



2 < C V / «^ +fe r 4 |d r fc u| 2 dr. (4.15) 



Similarly, if d k u G Y^ +fc /or eac/i < fc < f/3] + 3, then u G L°°(0, R). 

Proof. Since ||«||i<x>(o,ii) < CIMIi^ftVR) from the Sobolev embedding inequality, by combining 
with (|4.14p . we get the desired result. □ 

Lemma 14.41 and 14.51 can be directly used to obtain the embedding estimates for w p u for 
p > 0. For instance, if we apply Lemma 14.51 for wu in place of u and use (|4.1ip again, we get 
the corresponding embedding estimate for wu. Here we record the results. 

Lemma 4.6. Let p > be given. The followings hold. 

1. Let d k r u G X /3+p+fc for each < k < \/3] + 2 w/iere \(3] > 1. T/ien u> p / 2 u G L°°(0,R) 
with the following estimate: 



IK /2 «||;L < C V f w 0+p+k r 4 \d k u\ 2 dr. (4.16) 
fc=o ■ / ° 

Similarly, if d k u G y/ 3 +P+ fc for each < k < \/3]+3, then w p / 2 u G L°°(0, R). 
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2. Let d k u € X? +p+k for each < k < \/3] . Then w p / 2 r 2 u G L 2 (0, R), equivalently u £ X p 
with the following estimate: 



.8 m , R 

fe=0 

4.2 L°° bounds 



/ w p r A u 2 dr <CY^ w P+p+k r A \d k u\ 2 dr. (4.17) 



A direct consequence of the above Hardy embedding inequalities is the validity of the small- 
ness assumption (|4.7p within our energy space induced by £. 



Lemma 4.7. There exists C > so that 

r [<»l+3 j i ^ r r«i+ 3 i 

|C| + M+ £ h^C|+ £ |^fC,|<C^. 

q=l g=l 

Proof. We will present the detail on the following three terms 

Other terms can be treated in the same way. To see the boundedness of £ r , we apply Lemma 
3} take f3 = a + 1 in (|4.15p to deduce that 

N+3 



llCrllL <c"t I* w a+1+k Ad k r Q 2 dr < C£. 
k=0 Jo 



For w [L ^ ]/2 dt 2+ ]+1 C, we will apply Lemma 023 Take = a - [^^] in (13151) . Notice 
that when 3 < a < 4, f3 = a — 3; when 4<a<5, /3 = a — 4 and therefore, \0\ = 1, which 
satisfies the condition on f3 in Lemma 14.61 Then, we see that 



1/2 ^ 1+1 CllL < / w*+ k Ad k r d [ — ]+ \\ 2 dr 

k=o Jo 

<cJ2^^ ]+k ' k by m- 



k=0 



J t 

will apply Lemma 14.61 for the last term. Take /3 = a — [ — ] so that = 1 in ()4.16p . We 

\W [ 2 



Since [i^M] + 3 < [a] + 3 for a > 3, we conclude that H^ 1 ^/ 2 ^ 2 ^CIIL 3 £■ We 
will apply L 
first obtain 

rM+3i J rR r r a l+3, 

-^ 2 d [ — ] C r \\l <CV w a - 1+k r 4 \d k d [ — ] ( r \ 2 dr. 

k=0 J ° 

Now by applying the Hardy inequality (|4.11j) and by using the definition of the energy (14. 5p 
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we further deduce that 



4 fR rM+3, 

r 



) [^- 1 ]/ 2 ap 1 C r ||L<cf; I w a+1+k r^\d k d, 

k=0 J ° 



k=0 

Since [^M] + 4 < [a] + 3 for a > 3, we conclude that \\w [S ^- 1 V 2 df^ ] C r \\l 3 £• □ 
We are now ready to prove Proposition 14.11 the equivalence of £ (t) and £(t). 

4.3 Proof of Proposition 14.11 

Denote two energy terms in (|4.5p by 8( ,k and f^''* so that 

gj,k _ gj,k _|_ 1+a qt 

where 

£/> fc = [ R w a + k r A \d{- k d k ( t \ 2 dr and ^'' fc = /* u> 1+a+fc r 4 |^'- fc a*Cr| 2 dr. (4.18) 
Jo JO 

Notice that by definition, 

= £ j and £/' fc = £^ k ~ l for 1 < A; < j. (4.19) 

We will prove the equivalence for the simplest case: j = 1 and k = 1 first and then move 
onto the general case j > 2. 

Lemma 4.8 (£ ' )• There exists a constant C > suc/i i/iai 

f 1 ' 1 < C^+f 1 ). 

Proof. In this case, because of (|4.19j) . we only need to show that J ^ u> 2+c V 4 |C rr | 2 dr is bounded 
by the temporal instant energy. We recall the equation (j4. 1 j) in the following form: 

l{w 1+a r%) r = + (3 7 - 4)w a r^(r)C - w a r^(r)f(C) 

+ r 3 (w 1+a h(C,Crr)) r 

where /(C) = 0(|C| 2 ) and h(^,(^ r r) = 0(|C| 2 )+0(|Cr r | 2 )- For notational convenience, we have 
used 7 = (1 + a)/(x. We will exploit the elliptic structure of the term in the left-hand-side of 
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(|4.20p . Square both sides of (|4.20p . divide them by w a r and integrate it over (0, R) to get 



R -v 2 9 
—^—7 I (w 1+a r 4 ( r ) fdr< 
w r ^0 



H w a r 4| Ctt ,2 



dr+ / w a r 4 \$(r)f(()\ 2 dr 



(*) 



(0 



(«) 



r-R /■« i , 

+ / w Q r 4 \<S>(r)(\ 2 dr + j -^\r 3 (w 1+a h((,(rr)) r \dr 



(4.21) 



(iii) 



(it)) 



It is clear that 



(4.22) 



For (it;), we first see that 

r 3 (w 1+a h(C, ( r r)) r = w 1+a r A 8 2 h ( rr + w 1+a r 3 (d 1 h + Sb/OCr + r 3 (w 1+a ) T h 

where d\h and d 2 h mean the derivative of h with respective to the first and second argument 
respectively. For instance, if h(C,Cr r ) = C 2 + (Cr^) 2 , d\h = 2( and d 2 h = 2£ r r. By using the 
notation given (14. 2p . we write (tf 1+Q ) r = -r?o a $(r) and hence, we see that (iv) is bounded 
by 

pR pR pR 

(iv) ^ / w 2+a r 4 |a 2 /i| 2 |Crr| 2 *+ / ^ 2+Q r 2 |ai/i + a 2 /i| 2 |Cr| 2 ^+ / /r 4 1 $(r) 1 Vdr . 
Jo Jo Jo 



(tt))l 



(™)2 



It is easy to see that the first term in the right-hand-side is bounded by 9 2 £ r ' . For the second 
and third terms, we will employ Hardy inequalities; by using the Hardy inequality near the 
origin (|4.9j) for (iv)\ and the Hardy inequality near the boundary (I4.11|) for (iv) 2 and noting 
that |di/i + d 2 fr| 2 < 6 2 and \h\ < 9f, we deduce that 

(iv) < C9j{£° + S 1 ' 1 }. (4.23) 

We now turn our attention to the term (*) in the left-hand-side of (14.2ip . First notice that 

(w 1+a r 4 ( r ) r = w 1+a r 4 C„. + 4w 1+a r 3 ( r + (w 1+a ) r r 4 ( r 
= w 1+a r% r + 4w 1+a r 3 ( r - w a r 5 <5>(r)( r 

Thus the (*) in (14.2ip reads as 



pR 

(*) = 7 2 / w a r 4 
Jo 



wCrr + - r$(r)Cr 



dr. 
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By expanding terms out, we see that 

(*) rR rR rR 



7 2 



= f w 2+a r% rr \ 2 dr + l§ [ w 2+a r 2 \Cr\ 2 dr + [ 
Jo Jo Jo 

pR rR p 

+ 8 / w 2+a r 3 ( rr ( r dr-2 / w 1+a r 5 $(r)( rr ( r dr -8 / 
Jo Jo Jo 



w a r b \$(r)\ z \( r \ z dr 
Ft i*R pJR 



w 1+a r A $(r)\tr\ z dr 



(*)i (*)a 
For the (*)i and (*)2, we integrate by parts to get 

pR t-R 
= -4 / (w 2+a ) r r 3 \( r \ 2 dr - 12 / w 2+a r 2 \( r \ 2 dr 
Jo Jo 

fR 



9 4- n f H f 1 
A—— / w l+a r^{r)\Q r \ 2 dr- 12 / 
1 + " Jo Jo 

pR pR pR 

- w a r 6 \$(r)\ 2 \( r \ 2 + 5 / w l+a r^{r)\Cr\ 2 dr+ / u> 1+f V 5 $'(r)lCr| 2 dr 
Jo Jo Jo 



Hence we obtain 



w 2+a r 2 \( r \ 2 dr = { -4 

7 



pR pR. 
/ w 2+a r 4 \Q rr \ 2 dr + 4 / ■ 
Jo Jo 

/•« 2 , pR pR 

+ 3 / w 1+< V 4 $(r)|Cr| 2 dr-4^-^ / w 1+a r 4 $(r)\( r \ 2 dr - / w 1+a r 5 $'(r)lCr| 2 dr 
Jo 1 + « Jo Jo 

N v ' 

It is clear that the last three terms in the right-hand-side are bounded by the zeroth-order 
energy £°. Now combining it with (jOTj) . (j4T22|) . and P~23|) . we get 

For sufficiently small 9\ > 0, by absorbing the first term in the right-hand-side into the 
left-hand-side, we deduce that 

Sr' 1 <C 2 (£° + £ 1 ) 

for some constant C 2 > 0. This finishes the proof for the case of j = 1 and k = 1. □ 

We now turn into the cases j > 2. As in the case of j = 1, we will directly use the 
equation and take advantage of the elliptic estimates. What is subtle here is to capture the 
correct behavior of solutions in the normal direction d r near the boundary. 

Lemma 4.9 (£ J,fc ; 1 < k < j, 2 < j < \a\ + 3). There exists a constant C > such that 

j 

£ j,k <c^£\t). 

1=0 

Proof. Notice that because of (|4.19j) . it suffices to show that each spatial energy term £p k 
for 1 < k < j is bounded by £ (t). We will present the detail for j = 2; other cases follow 

2 1 

by the induction on j,k. When k = 1, the spatial energy term £ r ' contains one temporal 
derivative and two spatial derivatives. The time derivative of the equation (14.200 is a good 
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place to start. But then, the time derivative does not affect the weights at all since w and 
r do not change in time. Therefore, following the same procedure for E^ 1 in the previous 
lemma, we can deduce that 

8 2 > 1 <C{£ Q + E 1 + E 2 ) 



for some constant C > 
To deal with £ 2 ' 2 whi 

for Q rrr from (|4.20p . Here is the algorithm to do so: first divide both sides of (|4.20p by r 3 w a : 



2 2 

To deal with £ r ' which contains three spatial derivatives, we will first derive the equation 



-y{wr( rr + (1 + a)w r r( r + 4< r } = + ~ 4 )r$(r)C - r*(r)/(C) 

+ {wrd2h( rr + w{d\h + c?2^)Cr + (1 + a)w r h] 

Then we take d r of both sides of the above equation and move the terms involving £ r into 
the right-hand-side to get 

j{wr( rrr + (2 + a)w r r( rr + 5wCrr} = — t{(5 + a)w r Cr + (1 + a)w rr r( r ^ 

+ (TT|)2 + (TTcF " £t$S + & ~ 4)(r*(r)C) r " (r*(r)/(C)) r (4.24) 

+ (wrd2h( rr + w(dih + c^Kr + (1 + a)w r h) 

As in the previous lemma, we square both sides of (|4.24p . multiply by w 1+a r 2 - here the 
multiplier w 1+a has been chosen inspired by the analysis carried out in [8] - and integrate it 
over (0, R) to get 



/ w +ol r 2 \wr( rrr + (2 + a)w r r( rr + 5w( rr | dr 
Jo 

(*) 

R 

w 1+a r 2 (\w r \ 2 + \rw rr \ 2 )\Cr\ 2 dr 



-< 



w 

W (1 + C) 4 dr + /o ^TTO^ dr + y o (1 + C) 6 ^, (4.25) 

v s V ' 

+ / <u; 1+Q r 2 



(r$(r)C) r | 2 dr + / w 1+c V 2 | (r$(r)f(()) r \ 2 dr 
Jo 



(iii) 

-R 



r 

+ / w l+a r 2 | (w;r5 2 /iCrr + w{d x h + d 2 /i)Cr + (1 + a)w r ti) \ dr 
Jo 

V v ' 

(iv) 

Note that (i), (ii), (iii) contain stronger weights near the origin and for those terms, we can 
use the localized Hardy inequality (|4.9p to obtain 



<£° + £\ (ii^^+f 2 , (iii^^ + f 1 . (4.26) 
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The last term in the right-hand-side of (|4.25[) : (iv) includes (, rrr , Cm Cr with different weights 
and it can be treated in a similar way as done for (iv) of (|4.2ip in the previous lemma. We 
expand it out and apply the Hardy inequalities both near the origin (|4.9p and near the 
boundary (14.111) to deduce that 



(IV 



<e\{£ Q + £}> l + £ 2 ' 2 }. (4.27) 



What follows now is the elliptic estimate for the (*) in the left-hand-side of (|4.25p . which 



will give rise to the term £ 2 ' 2 . 



(*)= ^ w 3+a r A \( rrr \ 2 dr + (2 + a) 2 f* w 1+a r A \w r \ 2 \( rr \ 2 dr + 25 f* w 3+a r 2 \( rr \ 2 dr 
Jo Jo Jo 

j-R t-R 
+ 2(2 + a) / w 2+a r A w r ( rrr ( rr dr + 10 / w 3+a r 3 ( rrr ( rr dr 
Jo Jo 



(*)l (*)2 

-R 



10(2 + a) / w 2+a r 3 w r \C rr \ 2 dr 
Jo 



For the (*)i and (*)2, we integrate by parts to get 

<*R fR pR 



(*)l 
2 + a 



f {w 2+a ) r r A w r \Crr\ 2 dr -A [ w 2+a r 3 w r \C„\ 2 dr - [ w 2+a r A w rr \Q rr \ 2 dr 
Jo Jo Jo 

pR i-R 

(2 + a) / w 1+a r A \w r \ 2 \( rr \ 2 dr - 4 / 
Jo Jo 



w 2+a r 3 w r \(; rr \' z dr 



R 

w 2+a r 4 w rr | C rr \ 2 dr 

fR fR 
(*) 2 = -5(3 + a) / w 2+a w r r 3 \Crr\ 2 dr - 15 / w 3+a r 2 \C rr \ 2 dr 
Jo Jo 



Thus we obtain 

fR r R 



f u,3+« r 4|Crrr| 2 dT- + 10 / W 3+a r 2 \( rr \ 2 dr 

Jo Jo 

= (*) - (a - 3) / w 2+a w r r 3 \( rr \ 2 dr + (2 + a) / w; 2+a rVr|Crr|^r- 



By noting iu r = — r$(r)/(l + a), we see that the last two terms are bounded by £ r ' . By 
combining with P~25j) . (|4T26|) and (|4T27]) . we deduce that 

^ 2 ' 2 <C7(£: + £: 1 +£: 2 + £: r 1 ' 1 ) 

for some constant C > 0. Since £r ^ £° + S 1 by the previous lemma, the desired result 
follows and this finishes the proof of the case j = 2. Other cases can be done inductively: take 
d r derivatives of the equation (|4.24p . square it, multiply by appropriate weights depending 
on the number of spatial derivatives, and exploit the Hardy inequalities and the elliptic 
estimates. The procedure and the estimates are similar to the previous cases and we omit 
the details. □ 
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5 Nonlinear weighted energy estimates 

In this section, we develop the nonlinear energy estimates for sufficiently small perturbed 
solutions (Ct,C) to the Euler-Poisson system f|2. lTj) or (|4.ip . Because of a subtle nonlinear 
structure of the pressure gradient term, a splitting to a linear part and a nonlinear part as 
done in (|4.ip may destroy the important structure during the estimates and it may lose a 
necessary cancelation property unless one would be extremely cautious. Here in order to 
obtain the close energy inequalities, we will work with the original perturbed form (|2.17p 
and introduce the nonlinear high-order energy norms in (|5.7p . We first recast (12.17P in the 
following form: 



w a r i C,tt 
(1 + C) 2 



+ r 3 d r 



w 



l+a 



i + \(Ac + c 2 + c 4 



+ 



i-(i + C) 4 
(i + C) 4 



1 



w a r^{r) = 



(5.1) 



where we have used the Lane-Emden equation (|2.15p with the notation introduced in (14.2 
and the following identity: 



(1 + C) 2 (l + C + Crr) = (1 + 2C + C 2 )(l + C + Crr) 

c 3 . 



1 + 3C + Crr + 3C 2 + 2C( r r + ( 3 + ( 2 ( r r 



= l + ^(r 3 (C + C 2 + 

We are now ready to perform the energy estimates. Throughout the section, the smallness 
of the solution (|4.7p is assumed. We start with £ °. 

Lemma 5.1 (£ °). For any fixed small n > 



w a r\ 2 dr 



(5.2) 



where \K \ < CO^ . 

Proof. We begin by multiplying (|5.ip by (1 + C) 2 Ct an d integrating over (0, R): 



R 



a 4 

w r 



^«lr + I ^(1 + 2 (tr 3 d r (w l+a { [l + ^(r 3 (C + C 2 + £) 



}) 



dr 



+ 



[i-(i + 4 ]C* 



w a r 4 &(r)dr = 



(i + C) 2 



We denote the left-hand-side by (/) + (II) + (III). We will estimate it term by term. Note 
that the first term (I) forms a perfect time derivative: 



1 d 
2~dt 



{I) 

For (//), we first integrate by parts 

(//) = - j\r\l + C)\ t ) r w l+a { 



R 



w a r A \( t \ 2 dr. 



l + ^(r 3 (C + C 2 + ^' 



l+ct 



1 \dr 
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and use the identity (1 + () 2 Ct = (C + C 2 + to see that it also forms a perfect time 
derivative 



(II) 



R 



[AC + C 2 + C Y n 



)y +a { 



i + ^p(C + C 2 + ^ 



l + g 



} 



1 



(ft 



/* ^+«{ar 2 ([l + 1 (r 3 (C + C 2 + ft) ] ' ~ l) + (r 3 (C + C 2 + $ )) } 

«/0 7* x T 



dr 



By Taylor's theorem, we see that the inside of {•} is non-negative and it can be written as 



,1+a 



2 + f) = ^(l + -)|3rC + r 2 C r | 2 + f) 



a 



where f) together with the weight w +a is bounded by Q\S . Thus we may write 



2 a dt J Q 



w 1+a \3r( + r 2 ( r \ 2 dr -^K 
1 1 ere 



(*) 

where iT < Cfli£ . The term (*) gives rises to the other part of 8°. To see it, we expand it 
out: 



R 



w 



\3r( + r 2 ( r \ dr = 9 / w 1+a r 2 \(\ 2 dr + / w 1+a r 4 |Cr| 2 dr + 6 / w 1+a r 3 CCrdr 



For the last term in the right-hand-side, by integrating by parts we obtain 

pR pR pR 

6 / w 1+a r 3 (( r dr = -6 / w 1+a r 2 \Q\ 2 dr - 3 / (w 1+a ) r r 3 \C\ 2 dr 
Jo Jo Jo 



and in turn 

*R 



pR rR pR pR 

/ w l+a \3r( + r 2 ( r \ 2 dr= / w 1+a r 4 |Cr.| 2 dr + 3 / ^ 1+a r 2 |C| 2 rfr - 3 / ( 
Jo Jo Jo Jo 



{w 1+a ) r r 3 \(\ 2 dr 



Notice that all three terms are positive. 

Lastly, we use the Cauchy-Schwartz inequality to control (III)- Since <3?(r) is bounded 
and by using (|4.7j) . it is easy to deduce that 



pR pR 
\(III)\<k w a r%\ 2 dr + C K w a r\ 2 dr 
Jo Jo 



This establishes d). 



□ 



For higher-order energy estimates, we will look at the weighted quantity instead of dealing 
with Q and £ r directly: this way, we will have more effective way of getting the estimates. In 
fact, notice that the nonlinearity in (I5.ip is directly related to J defined in (|2.1ip : 



i + \(A( + t 2 + c 4 



(i + C) 2 (i + C + Crf) = £ 2 (£ + £, r r) = J 
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w a r 4 Cttt _ w"r 4 CaCt _ w a r 4 <f>(r)Ct 
(1 + C) 2 (1 + C) 3 (1 + C) 5 

- ^ r 3 d r (w 1+a [1 + ^(r 3 (C + C 2 + £))r] • ^ (r 3 (l + C) 2 C*) r ) = 



and hence, the equation (|5,ip can be written as 

- 1}) + L^+0! ^ (r) = 0. (5.3) 

Moreover, we observe that 

This motivates us to consider the following variable: 

tp=(l + C) 2 Ct so that J t = ^[r 3 <p] r ■ (5.4) 
The time derivative of (15.111 reads as 



(5.5) 



By using the definition of ip, we also obtain 

vt = (i + c) 2 c« + 2(1 + oc? and ^* = (i + o 2 c«t + 4(1 + octCtt + 2d 5 

and hence by a straightforward computation, one can see that (|5.5p reads in terms of p as 
follows: 

w°'r 4: (ptt 6w a r 4l (pipt 14w a r 4 p 3 4u> c Y 4 <3?(r)y? 

(ITO 1 " (1 + C) 7 " (1 + C) 10 " (1 + C) 7 (5 6) 

- l -±^r%. (w^J-^hr^} r ) =0. 
a \ r l ) 

The advantage of the above tp equation (|5.6j) is that the last term, which contains full spatial 
derivative, is self-adjoint and more or less linear with respect to (p up to lower order terms, 
and that such a structure will not be destroyed under the time differentiations. This makes 
the higher order temporal energy estimates affordable. 
We now introduce the nonlinear ^—energy <£ l for i > 1, 



e=l W V dr y dr + — f\ 1+a J-^\\(r s dt^)r\ 2 dr (5.7) 



H - f ^ 4 |arVt| 2 1 + a r R 

1 + C) 4 + a J " " 

When i > 1, (£* corresponds to the homogeneous weighted energy and we will show in 
Lemma 15.31 at the end of this section that they are equivalent under the assumption (|4.7|) . 
We record the high order energy inequalities for the solutions to (|5.6p . 



Lemma 5.2 (£ l ; 1 < i < \a] +3). Under the assumption (|4.7p . i/ie solutions to (|5.6p satisfy 
the following: for any small fixed k > 

1 d * 1 

< (C0i + + C£° + £ J . 

— ■ tit 
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Proof. We start with i = 1. We multiply (|5.6p by ^ and integrate it over (0, R). We denote 
each integral by Ik for 1 < k < 5. We will estimate them term by term. I\ forms an energy 
plus a commutator and thus 1% + 12 can be written as 



r r 1 d 

h+h =2dt 



R uj«7-4|y, t |2 

(i + C) 4 



(i + C) 5 



-<ir 



where we have used (I5.4p for 72- Note that the second term is bounded by O1& 1 since < Oi 
due to (|4.7p . It is easy to see that I/3I < C8\(<£ 1 +S°). For I4, <3?(r) is bounded and by using 
the Cauchy-Swartz inequality, we see that 

|/ 4 | < kC 1 + C k £°. 



We next move onto I5. We first integrate it by parts in r and then we see that it forms a 
perfect time derivative plus a commutator: 



h 



l + a 



R 







a 

11 + a d 
2 a di 



(r 3 ip t ) r ■ w 1+a J ^(r 3 ip) r dr 



in 



l+a j- 



' :m 1 1 , 3 N ,2 







1 1 + a 1 + 2a 

2 a a 



J" 



(r <^) r dr 



1 



Jt~? \(r <p)r\ dr 



Since J is bounded away from zero and from above and Jt is bounded by 6\ because of (|4.7p , 
the second term is bounded by OitB 1 . This completes the case of i = 1. 

For higher order estimates of (£* where 2 < i < \a] + 3, we first write out the (i — l) th 
temporal derivative of (|5.6p : 



rv 4 *~ti — 1 ^ 1 



(1 + C) 4 

ar 1 



i=i 
14u> Q rV n 



1 



3 



,1-1 



(i + C) 7 



(i + C) 



10 



i=i 



.1+0 



(i + C) 7 

l + 2a 1 



l^r 3 a r (V+ Q 1 [r 3 ar V] 



(5.8) 



r 3 ar 1_ V]r) = 



where cij and C2j are binomial coefficients. Now we multiply (|5.8p by <9j</? and integrate it 
over (0,-R). We denote each integral by for 1 < k < 7. As before, we will estimate them 
term by term. As in the case of ii, the first term J\ forms an energy plus a commutator that 
is bounded by 6\£ l : 



Jl ~ 2 dt 



(i + C) 4 



dr + 



R 



2w a r A Ct{d\(py 
(1 + C) 5 



dr 



Next, J2)^3)<^4 w ih yield lower-order nonlinear terms. Here we present the detail for J 2 . 
We may assume 1 < j < [|] + 1 since the other cases can be treated in the same way. We 
treat the case when dt hits ( j times for the second factor dfl ^^ ]. By Cauchy-Schwartz 
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inequality, we first get 

| H < Csuph^ U (1 ^V * U ) ' (5.9) 



(*)l 



«2 

Note that < 6\ by the assumption (|4.7jl . For (*)2, by Lemma [5T3], it suffices to consider 



R _ _ _ xl 



(*) 2 = / w a - J+1 r 4 |c>r J+2 C| 2 ^ 



o 



Then by the embedding inequality (|4.17p in Lemma 14,61 we get 

i- 1 r R J- 1 i- 1 



[(*) 2 ] 2 < I w a - j+l ™ 3 ~ l+k Adl- 3+2 d k r (;\ 2 dr < Y^£ l -i +1+k ' k < C^S 

k=0 fc=0 fc=0 

and hence, we deduce that 



\j*\ < c^i j^^-i+i+fc (r) 2 < ^ + 

\fc=0 / k=0 

The J5 is also lower-order. When dt hits ip (i — 1) times, we use the Cauchy-Schwartz 
inequality to derive that it is bounded by 

' ' " 1 -dr + C K / — ' * , ' dr 



o (l + c) 4 K 7 (i + C) 4 

and for other cases, we can use the above standard nonlinear estimates to conclude that 

i-2 



|J 5 | < (K+cei)e + c K e- 1 + ce 1 J2 sk 



k=0 

As in the previous I§ case, the Jq term will form an energy: 



a 







11 + a d f I , l + 2a 1 1 3„i_1 s |2 , 

' w i+a J - (r 5j V)r dr 



2 a di Jo 7,2 
1 1 + a 1 + 2a f R 1, i+3a 1 1 3 ,1 ,2 , 



2 a a 

It is clear that the second term is bounded by 9\<t % . It now remains to estimate J7. Like the 
J2 term, this is a lower order nonlinear term. We may assume 1 < j < [|] since the other 
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cases can be treated in the same way. By Cauchy-Schwartz inequality, we first get 



j — l 

\Jj\ < Csup \w~cP t J 
v 

(*)s 



(r 3 ar 1_j V)r 



dr 



(*) 4 



(5.10) 



Note that (*)3 < 6\ by the assumption (|4.7p . For (*)4, by Lemma [5.31 it suffices to consider 

1 



2+a -^r 4 \d z t - J Cr\ 2 dr 



Then by the embedding inequality (14. 17ft in Lemma I4.6 



[(*) 4 ] 2 < E 



k=0 



w 2+a ' j w j+k r 4 



fc=0 



A;=0 



and hence, we deduce that 

i 

/ 3 \ 2 i 

|J 7 | < C0i < Ox? + c9iJ2s i - j+k . 



\k=0 



k=0 



□ 



Next we establish the equivalence of energies (5: and E under the assumption ()4.7f) . 



Lemma 5.3. Assume (|4.7j) . For eac/i 1 < i < \a\ + 3 ; there exists a constant C > smc/j 

(5.11) 



e = e i + # 



w/iere |^| < C0i Efc=o f ■ 

Proof. This can be established by direct computation with the aid of Hardy inequality (|4,9 
and Lemma 14.91 We will prove the case of i = 1. Prom (|5.4p and (|5.7p . 

rB ^ r 4| (1 + C) 2 Ctt + 2(1 + c)c 2|2 



!i + C) 4 



-dr 



+ ^/ ^j-^4l(- 3 d+c) 2 c t ),| 2 * 

a Jo r 
= (') + («) 
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Now 



(/)= / m <- r >\( u + ^L\ 2 dr 

\„V M V + [ + [ f£t ir S El + M 

and 

TT-W= [ R w 1+a r 4 J- 1 ^(l + 0*\dt(r\ 2 dr+ r^J-^l|(r 3 (l + C) 2 Ut| 2 ^ 
1 + OL Jo Jo r 

+ 2 f R w 1+a J-^r(l + C) 2 0r(r 3 (l + C) 2 U^ 

JO 

so that 

4 = J R w 1+a r' [ J- 1 ^ (1 + C) 4 " 1] \dtCA 2 dr + jf * w 1+a 1 | (r 3 (l + C) 2 U* |* dr 

+ 2 J-^r(l + C) 2 C*,(r 3 (l + CfUtdr 

Jo 

Hence, we can write (S 1 = (I) + (II) = S 1 + R\ + ^-^^\- It is easy to see that 

For we denote three integrals by (i) + (ii) + (in). Notice that 

I (i) I 2*1 [ R w 1+a r 4 \d t C r \ 2 dr since | J~^(l + C) 4 - l| 3 9 X 
Jo 

l(*OI 3 *l (j£ w 1+a r 4 \(r\ 2 dr + u» 1+Q r 2 |C| 2 ^ since |Ct| < 9\ 



< ft? 



ft\ {^j w 1+a r 4 \C r \ 2 dr + J w 1+a r 4 \C\ 2 dr^j by Hardy inequality ([M 

^ 

1+«J|A |2j„ , / „„l+ai\/- i2 . , / „.,l+a„4|A|2, 



|(m)|^0i / w i+Q r 4 |Ct r rdr+ / u/ i+a r 4 |Cr-r<ir + / u; i+c V|Crdr 
VJo Jo Jo 

where in order to obtain the estimate of (to), we have applied Cauchy-Swartz inequality and 
used \Ct\ < 0\ and the estimate of (ii). Therefore, we conclude that 

iMl^c^ + f ) 

and this finishes the proof of i = 1. 

Other cases of i > 2 can be treated in the same manner by expanding d t t ^ 1 ip in terms of 
C, £t,...,c^C in the energy form (|5.T|) : since 

i-1 

^- v = (i + cfdic + £ c ^'[d + o 2 ] a^c 

i=i 
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where c,'s are binomial coefficients. By plugging into (|5.7p . we see that the leading-order 
term gives rise to £ l . Following the same spirit of the above computations on &\ and ^ 
and exploiting the nonlinear estimates as done in the previous lemma by using (|4.7|) . Hardy 
inequality (I4.9P and Lemma l4.9l one can deduce that lower-order nonlinear terms are bounded 



6 Nonlinear instability 

Based on the nonlinear estimates in the previous section, we are now ready to show a boot- 
strap argument that allows us to control the growth of £ (t) in terms of the linear growth rate 
y/JIo, showed in Sectional The idea is to assume that the lowest order energy £° starting from 
small initial data grows no faster than the linear growth rate; then the energy inequalities in 
the last section allow for a bootstrap argument that claims that all of £ grows no faster than 
the linear growth rate. 

Proposition 6.1. Let (dtC,, C) be a solution to the Euler-Poisson equation (|4.ip such that 
y/£(0) < C d and >/f°(t) < CoSe^ 1 for < t < T 

where £°(t) and £{t) are as defined in (|4.3p and (|4.4p . Then there exist C+ > and 6+ > 
such that ifO<t< min{T, T(S, 0*)} then 



where we have written T(5, 9+) = 




Proof. We will employ a bootstrap argument using all of the nonlinear energy estimates 
derived in Section [SJ We now choose 0\ and k sufficiently small (but fixed) in all of the 
estimates in Lemma [5.21 and [5.31 so that C9\ + k < and CB\ < 1/8 in all of the energy 

inequalities. As before, throughout this proof we will write C for a generic constant; we write 
this in place of C to distinguish the constants from those appearing in the nonlinear energy 
estimates. 

Define T* by 

T* = sup{i : V£(s) < 0i for < s < t} 

Then by Lemma 14.71 the assumption (|4.7p is satisfied by the solution to the Euler-Poisson 
equation ([43]) for < t < min{T,T*}. 

Let < t < min{T, T*} be given. From Lemma [5.2l we have that for each 1 < i < \a \ +3 

±d i < yW^ C£0 c y j ( } 

2dt ~ 2 ^ y ' 

We begin with i = 1. Notice that for i = 1, the last term in ([6.1|) is not present. Since 
£°(t) < CQ5 2 e 2 ^ I ° t by the assumption, by using Gronwall inequality, we see that 

< C5 2 e 2 ^ t . 

By Lemma EH we deduce that 

£ l < C8 2 e 2 ^ t . 
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For 2 < i < \a] + 3, inductively, we apply Gronwall inequality to (|6.ip to obtain 

e < C5 2 e 2 ^ )t 
and in turn by using Lemma 15.31 we derive that 

£ i < C5 2 e 2 ^ t 

for all 1 < i < |~a] +3. Moreover, by Lemma [4, 9 \ we also deduce that for < t < min{T, T*}, 

£{t) < C5 2 e 2 ^ t (6.2) 

for some constant C > independent of 5. We now choose 0* > such that C(0*) 2 < (0i) 2 . 
We consider the following two cases. 

(i) T(5,9±) < min{T, T*}: In this case, the conclusion follows without any additional 
work. 

(ii) T(6,6+) > min{T,T*}: We claim that it must hold that T < T* < T(S,6 it ) in which 
case the conclusion directly follows. To prove the claim, we note that otherwise we would 
have T* <T < T(5, 0*). Letting t = T* from ([62]), we get 

£(T*) < C5 2 e 2 ^ T * < C5 2 e 2 ^ T ^ = C(9*) 2 

by the definition of T(5, 9±). However, this is impossible due to our choice of 8+ since it would 
then contradict the definition of T* . Since we then find our desired estimate in both cases, 
this concludes the proof of the proposition. □ 

In order to establish the nonlinear instability, we will take the growing mode profile (j)Q 
constructed in Section [3] as the initial data to the Euler-Poisson equation as in (j6.5[) . By 
Lemma |3.6| we deduce that all the spatial energies in £ for <j)Q are finite. See also Remark 
13.71 Then, of course, the initial temporal total energy £(0) for such <pQ is obtained from the 
initial spatial energy through the equation. 

We are now ready to prove our main result. 

Theorem 6.2. There exist 9q > 0, C > 0, and < 5q < 29q such that for any < 5 < 5q, 
there exists a family of solutions (dtC S , C 5 ) to the Euler-Poisson system (|4.ip so that 

y/£(0) < C5 but sup V£°(t) > O . 

0<t<T 6 

Here T s is given by T s = In 2 ^ L . 
Proof. We write the equation (14. ip as 

w a r% t = L( + N( (6.3) 
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where 

L( = l±^( w 1+a r\ r ) r + - 4)(1 + a)w a r 3 w r ( 

a a 

N( = (2C + C 2 ) l^^(w 1+a r%) r + (3-^ - 4)(1 + a)w a r 3 w r ( 
[a a 

- (1 + C) 2 r 3 d r (w 1+a HC, Crr)) - (1 + C)4 ~ 1 ~ 4 2 C(1 + C)4 r 3 a,(^ 1+ ") 

(l + Cr 

where /i is given by ()3.ip . Recall that L£ was introduced in (j3.5j) . Let 4>o be a growing mode 
constructed in Proposition 13.11 Normalize (fro such that 

(1 + tio)Ho\\x- + UoWy- = 1 ( 6 - 4 ) 

where || • and || • ||ya are given in (|3.19p . Consider a family of initial data ((,(t)\ t _ = 
^(^0) \fjM)4>())- Then due to the regularity of the growing mode 4>q obtain in Lemma [331 and 
\EM V£(0) < C 5 for some constant Co > 0. Let (C S Xt) be £— solutions to the Euler-Poisson 
equation (|6.3p with the given initial data 

(C S Xt)\ t=0 = S(4>o,^o)- (6.5) 

Define T by 

T = sup{r : £°(t) < 35 2 e 2 ^> t for < t < r}. 
Then by Proposition [67TJ there exist C* > and 9* > such that for < t < min{T, T(<5, 6*)}, 

y/£(t) < cje^ 1 < c*0* (6.6) 

We now wish to write the solutions as linear and nonlinear parts. To do so, let C{t) (^) 
be the solution operator for the linearized equation 

Ctt - ~^L( = 

with the given data (£)| t=0 = (^). Note tnat is S iven b ^ in (HUH])- Then 

by Duhamel's principle, the solutions to the nonlinear Euler-Poisson equation (|6,3p can be 
written as 

For notational convenience, we introduce the following norm: 

H*llo = ll^llxc + + II*IIk« 

We remark that || • ||q corresponds to the zeroth-order energy £°. Now since cj)Q is a largest 
growing mode, we see that 
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Then with (16.41). we obtain that 



C(t) 



Wo \ 



Serf**. 



(6.7) 



For the nonlinear part, we use Lemma 13.81 and 13.91 to derive that 



I C{t-s)(^_ 

JO \ w a r 







ds 



o Jo 



1 



:NC 5 (s) 



ds 



for a constant C > only depending on fi . Notice that N( = 0(\(\ 2 + \( r \ 2 + \( rr \ 2 ). We 
will show that || i 4 iV^(s)||x» is bounded by {5erf^ s ) 2 . To do so, first we crudely expand 
1 , NC S - 



1 



-NC 



(2C + C 2 ) 



1 + a(wl+Qr y- + (3±±^ - 4)(1 + a)^C 



n 



:i + o 



2 ^.(«; i+a ^(c, c,r)) (1 + o 4 - 1 - 4c(i + cr 



w a r 



a+o 5 



(1 + a) — 

r 



and then use the triangle inequality to get 
1 



-NC 



< C(|C| + IC| 2 )(IKrr| + \WrCr\ + \~^\ + |"^|) + C(| — + -^|). 



Then by using the L°° estimate in Lemma 14.71 and by the property of w in Lemma 13. 3\ we 
obtain 



1 



x a Jo 



R 



a 4 

to r 



^-T^C 6 (S) 
w a r 4 •> V ) 

R 



dr 



<C£(s)\ [ R w a+2 r A \C S rr \ 2 dr + f w a r 4 \(*\ 2 dr + [* ' w a+2 r 2 \C, 5 r \ 2 dr 
■/ io io 



(ai) 



(bi) 



/■fl fR fR . 

+ / w a r 4 \( s \ 2 dr + / to a r 4 |/t 5 | 2 dr + / w^Vl^/i^dr f. 
JO Jo Jo > 
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The first terms in the second and third lines have the right forms of the total energy £. For 
other terms, we apply the Hardy inequality near the boundary (|4.1ip to (a\) and (02) and 
the Hardy inequality near the origin (14.91) to (bi) and (62) to deduce that 



-^—[N( 5 (s) 

w a r 4 •> V J 



< C(£(s)Y 



x>- 
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Now by (|6.6p . we get 



C(t - s) 



ds 



<C I e 
o Jo 



£(s)ds 



< Cb 1 I e ^ (t - s) e 2 ^ s ds 



fl U2 



(6.8) 



< Ci{5e 



where C\ is a constant independent of 5. Now if necessary, fix 8* sufficiently small so that 
C X B* < 1/2. 

Next, we claim that T(5,8+) < T. Suppose not: T(5,9+) > T. Then by (|6.7p and (|6.8p . 
we have 



C(t) 



6yfjl(f)0 



(T) + 



f C{t-s) 
Jo 



ds 



(T) 



which would contradict the definition of T. With T(<5, 8+) < T in hand, we now see that by 
using (|6.7p and (|6.8p again, 

£0(T((5,e*)) = ||C*|| o Cr(5,0*)) > 5e^ TiS ' 9 ^ - Ue^ T ^ = -8*. 



Set 8q = ThO* and T s = T(5, 6+) = T(5, 28q). This finishes the proof of the theorem. 



□ 
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